Hyperon stars at finite temperature in the Brueckner theory 
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We perform Brueckner-Hartree-Fock calculations of hypemuclear matter at finite temperature and provide 
convenient analytical parametrizations of the results. We study then the properties of (proto)neutron stars con- 
taining hyperons. We find important effects of trapping and finite temperature on the structure of hyperonic 
stars. 
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I. INTRODUCTION 

The successful simulation of supernova explosions and the 
subsequent protoneutron star (PNS) evolution is still an open 
problem and currently a lot of theoretical activity is dedicated 
to it 1 1-9]. The fundamental input to these calculations is the 
nuclear equation of state (EOS) over a wide range of densi- 
ties, apart from microscopic information regarding diffusion 
and cooling processes. The output are time-dependent radial 
profiles of the thermodynamic quantities of interest, such as 
temperature, entropy, particle fractions, etc.. 

We have in previous articles investigated the nuclear EOS 
at zero and finite temperature within the Brueckner-Hartree- 
Fock (BHF) theory, which is currently one of the most ad- 
vanced microscopic approaches to the EOS of nuclear matter 
lfTol - [T4ll . We used the finite-temperature BHF EOS to model 
PNSs in our previous papers [ 15l-| l 8tl . limiting ourselves to in- 
clude hyperons as a Fermi gas ||15|]. In this work we further 
extend our approach by including consistently interacting hy- 
perons at finite temperature, and explore the consequences for 
PNS structure. We present in this article our results for the 
structure of cold NSs and hot PNSs containing hyperons. 

This is a technically demanding problem, which requires 
substantial numerical effort. Also due to this reason, hot hy- 
pemuclear matter has so far not been widely studied in the 
literature. We mention, however, the pioneering work regard- 
ing PNS evolution of Ref. |3], using a relativistic mean field 
(RMF) model including hyperons; a first investigation within 
the BHF approach [19], and a recent tabulation of hypemu- 
clear matter properties at finite temperature within the RMF 
approach |20ll. extending the finite temperature nuclear EOS 
of Ref. 

The PNS represents the typical state of the stellar object for 
some tens of seconds after supernova collapse, during which 
the system first deleptonizes and heats up the interior parts 
of the star in the process, before beginning to cool down by 
further neutrino diffusion. We do not intend to perform dy- 
namical simulations, but focus on the consistent construction 
of the temperature-dependent nuclear EOS and the evaluation 
of its basic consequences during the prominent PNS stage. 
We therefore assume strongly idealized, static profiles repre- 
senting this environment, namely we use a constant entropy 
per baryon throughout the star and investigate the sensitivity 



of the results to the chosen value of entropy SjA, as is often 

done iin Sin mm. 

We provide a short overview of the theoretical framework in 
Sect. II, before presenting our results in Sect. Ill, and drawing 
conclusions in Sect. IV. 



II. FORMALISM 

A. Brueckner theory at finite temperature 

The central quantity in the BHF formalism is the G-matrix, 
which in the finite-temperature extension lfTollT4ll24ll is deter- 
mined by solving numerically the Bethe-Goldstone equation, 
written in operatorial form as 
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(1) 

where the indices a,b,c indicate pairs of baryons and the Pauli 
operator Q and energy E determine the propagation of in- 
termediate baryon pairs. In a given baryon-baryon channel 
c — (12) one has 



6(12) = [l-ni(fci)][l -"2(^2)], 



(2) 



£■(12) =mi+t7i2 + ei{ki)+e2{k2) (3) 

with the single-particle (s.p.) energy e, 
the Fermi distribution = (gl^-C-O-A]/^ + l) the start- 
ing energy W, and the two-body interaction (bare potential) 
V as fundamental input. The various single-particle potentials 
within the continuous choice are given by 

t/l(^l)=Re E"(^2)(fcl^2|G(i2)(12) [£(12)] |^1^2)^, 

(4) 

where ki generally denote momentum and spin. For given par- 
tial densities p,- (/ = n,p,A,T.) and temperature T, Eqs. ([T]|4|i 
have to be solved self-consistently along with the equations 
for the auxiliary chemical potentials /I,-, 



Pi = Ln/W ■ 

k 



(5) 
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Regarding the interactions, in our calculations we use the 
Argonne Vig nucleon-nucleon potential |25] together with the 
phenomenological Urbana nuclear three-body forces (TBF) 
[26] . The corresponding zero-temperature nuclear EOS repro- 
duces the nuclear matter saturation point correctly and fulfills 
several requirements from the nuclear phenomenology ifTTi - 
[13I1 . In the hyperonic sector we employ the Nijmegen soft- 
core NY potentials NSC89 IH] fitted to the available experi- 
mental NY scattering data, see Refs. 1281 - 130(1 for details of the 
zero-temperature calculations. It turns out that at zero tem- 
perature only A and hyperons appear in the neutron star 
matter up to very large densities. We therefore restrict also 
the present study to these two hyperon species, neglecting the 
appearance of thermal iP and E+. 

Once the different s.p. potentials for the species / = 
A,E^ are known, the free energy density of hypernuclear 
matter has the following simplified expression 



TABLE I: Fit parameters for the free energy density, Eqs. | |12I19I >. 



(6) 



where 



Si = (ni{k)\nni{k) + [1 - ln[l - «,(^)]) (7) 

A- 

is the entropy density for component / treated as a free gas 
with s.p. spectrum ei{k) iflolfTill . 

A further simplification can be achieved by disregarding the 
effects of finite temperature on the single-particle potentials 
Ui, and using the T = Q results in order to speed up the cal- 
culations (frozen correlations approximation). This was the 
procedure followed in our previous publications fHKT^, and 
we apply it also in this work, due to the large number of cal- 
culations necessary when including the hyperonic degrees of 
freedom. 

All thermodynamic quantities of interest can then be com- 
puted from the free energy density, Eq. ^\ namely, the "true" 
chemical potentials Hi, pressure p, entropy density s, and in- 
ternal energy density e read as 

<3p, ' 
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(8) 
(9) 
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where p = p, is the baryon number density. We stress that 
this procedure fulfills by construction the Hugenholtz-Van 
Hove theorem in the calculation of thermodynamical quan- 
tities. For an extensive discussion of this topic, the reader is 
referred to Refs. ifTolfTill . and references therein. 



B. Parametrization of the free energy density 

The large number of degrees of freedom (4 partial densi- 
ties + temperature) renders inconvenient the use of the result- 
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ing hypernuclear EOS in tabular form. We therefore tried to 
approximate the numerical results by a sufficiently accurate 
analytical parametrization. We find that the following func- 
tional form provides an excellent parametrization of the nu- 
merical data for the free energy density in the required ranges 
of nucleon density (0.1 fm^^ ^ Pw ^ 0.8 fm^^), hyperon frac- 
tions (0 < Pa/Pn < 0.9, < Pz/Pn < 0.5), and temperature 
(0MeV<r<50MeV): 



f{p„,Pp,PA,Pz,T) =FnPn 
+ {Fa + Faa + Fal)Pa + 

+ {Fz + Fzz + Fza)PZ' 
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with the parametrizations at zero temperature: 

Fn = (1 -j3) (floPiV + ^oP^'') +i3 (fliPiV + ^iPj^' 



(12) 



(13) 



Fy = {a'^Y+a^x + alx^)pN + {b'^Y+l^YX + bYX^)PN ' (^4) 



Fyy' = aYY'p'iJ'" PyV , 
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(15) 
(16) 



where Pn = Pn + Pp\x = Pp/pN\ P = {l-2x)^; Y,Y' = A,E, 
and C = (3/5) (3;r^)^/^ « 5.742. At finite temperature the ex- 
pressions are extended as follows: 

Fn = Fm{T=0) 

+ [aot^pN + {dot^ + eoP)ln{pN)+fQt^/pN] (1 -j3) 

+ [ait^pN + {dit^ + eiP)ln{pN)+fit^/pN]P , (17) 
Fy = Fy{T^O) 

+ {dYt'- + eYt^)ln{pN)+fYt^/pN + 8Yt^HpY), (18) 
My = MyIT =Q)+bYt^pfj , (19) 

where t = r/(100MeV) and / and p, are given in MeV fm^^ 
and fm^^, respectively, (and iuax in MeV^'fm^^). 

Technically, these parametrizations were obtained by per- 
forming about 10^ BHF calculations at zero temperature in the 
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FIG. 1 : Free energy per baryon, F /A, at fixed nucleon density p/y = 0.6 fm^"^ and lambda fraction Pa/ Pn = 0.3, as a function of proton fraction 
Pp/Pn = 0, ...,0.5 and sigma fraction p^/pN = 0,0.1,0.2,0.3,0.5 for different temperatures T = 0, 10, ...,50 MeV. BHF data (symbols) and 
fit (curves) are shown. 



(p„,Pp,PA,P£)-space and then using the frozen correlations 
approximation to generate finite-temperature results, increas- 
ing by about one order of magnitude the number of "data" 
points /(p„,Pp,pA,pi,r). The fit parameters were then de- 
termined hierarchically for cold nuclear matter, cold hypernu- 
clear matter, hot nuclear matter, hot hypernuclear matter, so 
that the fits are optimized also in the more constrained cases. 
The final overall r.m.s. deviation of fit and BHF data points 
for F jA— f I p v& less than 2 MeV, which we consider fully 
satisfactory for our cuiTent purposes. 

As an illustration, we display some representative results 
for FjA in Fig. [l] namely a comparison of BHF data (sym- 
bols) and fit (curves) for fixed nucleon density = 0.6 fm^^ 
and lambda fraction Pa/Pn = 0.3, while varying proton frac- 
tion Pp/Pn = 0, . . . ,0.5, sigma fraction Pz/Pn = 0, . . . ,0.5, 
and temperature T = 0, . . . , 50 MeV. These are typical relevant 
values sampled in the parameter space of beta-stable hypernu- 
clear matter, as shown below. We notice an overall increase 
of the free energy with increasing fraction, for fixed T, 
which is due to the repulsive character of the effective T.^N 
interaction at this density. 



C. EOS of hot stellar matter and (P)NS structure 

In neutrino-trapped beta-stable nuclear matter, the chemical 
potential of any particle i = n,p,A,'L,l is uniquely determined 
by the conserved quantities baryon number B,, electric charge 
Qi, and weak charges (lepton numbers) ilf^ , : 



Hi = B,pLn - Qiilln -llp)+ if Vv, + 4'' Vv„ • (20) 



At given baryon density p = Y^tBiPi, these equations have to 
be solved together with the charge neutrality condition 



i 

and those expressing conservation of lepton numbers 



l^e,ll. 



(21) 



(22) 



As in our recent work lllTll . we fix the lepton fractions to Ye = 
0.4 and Y^^—Q for neutrino-trapped matter and treat the van- 
ishing of trapping in low-density matter ("neutrino sphere") 
llH I22I1 in an approximate manner. As in that reference, at 
subnuclear density, p < 0.1 fm^^, our BHF EOS is joined 
with the low-density finite-temperature EOS of Ref. |21] that 
accounts for clusterization of the matter, where the BHF ap- 
proach breaks down. 

The baryon chemical potentials required in Eq. (l20b are ob- 
tained from the free energy density /, Eq. (|6]l, and the chem- 
ical potentials of the non-interacting leptons from the free 
Fermi gas model at finite temperature. From the composition 
of beta-stable stellar matter, one can compute the total pres- 
sure p — ps + pi^, Eq. (|9ll, and the EOS p{e), with e = Cb + 
the total internal energy density, Eq. (fTTI) . The stable configu- 
rations of a (P)NS can then be obtained from the well-known 
hydrostatic equilibrium equations of Tolman, Oppenheimer, 
and Volkov [31] for pressure p{r) and enclosed mass m{r), 



dp 
dr 
dm 
"57 



Gme ( 1 + p/e) ( 1 + Anr'p/m) 



1 — 2Gm/ r 



(23) 



(24) 



(G is the gravitational constant). For a given central value of 
the energy density, the numerical integration of Eqs. ( |23l l and 
( l24l i provides the mass-radius relation. 
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FIG. 2: Relative particle fractions as functions of baryon density 
in beta-equilibrated matter at entropies S/A = 0, 1,2 without (upper 
panels) and with (lower panels) hyperons. 
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FIG. 3: Same as Fig.|2l but for neutrino-trapped matter. 
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FIG. 4: Pressure as a function of energy density for beta-equilibrated 
cold matter (solid curves) and neutrino-trapped hot matter at en- 
tropies 5/A = 1,2 (broken curves), without (upper curves) and with 
(lower curves) hyperons. The case with free hyperons at T = (green 
curve) is also displayed. 



TABLE II: Characteristics of the maximum mass configurations for 
different stellar composition and entropy. 
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Figs. |2] and [5] show the particle fractions at entropies S/A ~ 
0,1,2 in untrapped and trapped matter, respectively. We ob- 
serve the following qualitative features: (i) Finite temperature 
removes any particle thresholds, i.e., hyperons and leptons 
become more and more abundant at low densities with ris- 
ing temperature/entropy, (ii) Hyperon fractions are lower in 
trapped matter than in untrapped matter, in particular the 
is strongly suppressed, because due to the trapping condition 
it cannot easily replace the electron as is the case in untrapped 
matter, (iii) Compared to our previous work employing non- 



interacting hyperons (tist), one notes a slightly earlier onset 
and higher concentrations of the A, whereas the is a bit 
stronger suppressed. These properties are due to the attrac- 
tive/repulsive character of the effective A/E^-nucleon interac- 
tion in dense matter obtained with the NSC89 potential |29]. 

These features have direct consequences for the EOS p{£) 
that is shown in Fig.|4]for different configurations represent- 
ing cold untrapped NS matter and hot trapped PNS matter: 
For purely nucleonic matter the effects of trapping and tem- 
perature are not very large, but both soften the EOS. On the 
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FIG. 5: Gravitational mass (in units of tiie solar mass Mq = 1.98 x 10 g) as a function of radius (left panel) and central baryon density 
(right panel) for cold NSs (solid curves) and neutrino-trapped PNSs at entropies S/A = 1,2 (broken curves). The thick (thin) curves describe 
configurations with (without) hyperons. 



contrary, hyperons soften the EOS of untrapped matter much 
more than that of trapped matter, due to their higher concen- 
tration in the former. Altogether, finite entropy and in partic- 
ular trapping affect hyperon-rich matter much more (and in 
an opposite sense) than nuclear matter. For comparison, we 
display in Fig. |4] also the EOS for untrapped matter at T — 
obtained with non-interacting hyperons (green curve), which 
turns out to be very soft. We remind the reader that such an 
EOS gives a very low value of the NS maximum mass around 
IMollil. 

The relation p{£) as input to Eq. ( |23] ) determines directly 
the mass - radius relations of (P)NSs shown in Fig. |5] (left 
panel). Consistent with Fig. |4] one observes for nucleonic 
stars a slight reduction of the maximum mass (from about 
1.84 M0 to 1.74 M0) due to trapping and finite temperature, 
while for hyperon stars both trapping and also finite tempera- 
ture increase notably the maximum mass (from about 1.31 M0 
to 1.58 M0). The latter feature would permit a delayed 
collapse phenomenon, as is usually found for hyperon stars 
iSSIsllSJl- However, this conclusion is rather academic, be- 
cause the maximum mass of hyperon stars is 1.31 Mq in our 
approach, so that most observed NSs [33] would actually be 
hybrid stars involving a transition to quark matter in the inte- 
rior, as has been investigated in ll34ll . This is also pinpointed 
by the mass - central density relations, displayed in Fig. |5] 
(right panel), which shows for hyperon stars central densities 
up to about ten times normal nuclear matter density, where 
a realistic description of stellar matter should necessarily in- 
clude quark matter degrees of freedom. 

Table summarizes our results for the maximum masses 
of the different stellar configurations. As far as the minimum 
mass of PNSs is concerned, we find values slightly above 
0.5 Mq, thus confirming our results of Ref. 1 17], with a small 



discrepancy for the S/A — 2 case, which is due to the use of 
the frozen correlations approximation in the present calcula- 
tions. 



IV. CONCLUSIONS 

Summarizing, we have presented a convenient parametriza- 
tion of the free energy density of hypernuclear matter at finite 
temperature obtained consistently within the BHF framework 
using the VigH-UIX nucleon-nucleon and the NSC89 nucleon- 
hyperon interactions. 

Applied to the computation of (P)NS structure with sim- 
plified temperature profiles, we obtain relatively large effects 
of trapping and finite temperature in hyperon stars. However, 
their maximum mass is quite low, implying the presence of 
quark matter in the interior of heavier objects. 

For the future we hope to use improved NY potentials as 
well as TBF within the presented formalism in order to verify 
this important conclusion. 
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